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ABSTRACT 

A model of supernovae feedback during disc galaxy formation is developed. The model 
incorporates infall of cooling gas from a halo and outflow of hot gas from a multiphase 
interstellar medium (ISM). The star formation rate is determined by balancing the 
energy dissipated in collisions between cold gas clouds with that supplied by super- 
novae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created 
by thermal evaporation of cold gas clouds in supernovae remnants, and criteria are 
derived to estimate the characteristic temperature and density of the hot component 
and hence the net mass outflow rate. A number of refinements of the model are in- 
vestigated, including a simple model of a galactic fountain, the response of the cold 
component to the pressure of the hot gas, pressure induced star formation and chem- 
ical evolution. The main conclusion of this paper is that low rates of star formation 
can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with 
circular speed 50 kms~^ can expel ~ 60-80% of its gas over a time-scale of ~ 1 
Gyr, with a star formation rate that never exceeds ~ O.lMo/year. Effective feedback 
can therefore take place in a quiescent mode and does not require strong bursts of 
star formation. Even a large galaxy, such as the Milky Way, might have lost as much 
as 20% of its mass in a supernovae driven wind. The models developed here suggest 
that dwarf galaxies at high redshifts will have low average star formation rates and 
may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous 
discs might explain the numbers, metallicities and metallicity dispersions of damped 
Lyman alpha systems. 



1 INTRODUCTION 



Since the pioneering paper of White and Rees (1978), it has 
been clear that some type of feedback mechanism is required 
to explain the shape of the galaxy luminosity function in hi- 
erarchical clustering theories. The reason for this is easy to 
understand; if the power spectrum of mass fluctuations is ap- 
proximated as a power law P{k) oc k", the Press-Schechter 
(1974) theory for the distribution of virialized haloes pre- 
dicts a power law dependence at low masses 



dN{m) 
dm 



-(9-n)/6 



(1) 



For any reasonable value of the index n (n ~ — 2 on the scales 
relevant to galaxy formation in cold dark matter (CDM) 
models), equation (0) predicts a much steeper mass spec- 
trum than the observed faint end slope of the galaxy lumi- 
nosity function, dN{L)/dL oc L", with a ~ — 1 (Efstathiou, 
Ellis and Peterson 1988, Loveday et al. 1992, Zucca et al. 
1997). Furthermore, the cooling times of coUisionally ionised 
gas clouds forming at high redshift are short compared to 
the Hubble time (Rees and Ostriker 1977, White and Rees 
1978). Thus, in the absence of feedback, one would expect 
that a large fraction of the baryons would have collapsed at 
high redshift into low mass dark matter haloes, in contra- 
diction with observations. 



In reality, there are a number of complex physical mech- 
anisms that can influence galaxy formation and these need 
to be understood if we are to construct a realistic model of 
galaxy formation. In the 'standard' cold dark matter model 
(i.e. nearly scale invariant adiabatic perturbations), the first 
generation of collapsed objects will form in haloes with 
low virial temperatures (T ^ lO'^K, characteristic circular 
speeds Vc ^ 20kms^^). Molecular hydrogen is the dominant 
coolant at such low temperatures and so an analysis of the 
formation of the first stellar objects requires an understand- 
ing of the molecular hydrogen abundance and how this is in- 
fluenced by the ambient ultraviolet radiation field (Haiman, 
Rees and Loeb 1997, Haiman, Abel and Rees 1999). As the 
background UV flux rises, the temperature of the intergalac- 
tic medium will rise to ~ 10'*K {e.g. Gnedin and Ostriker 
1997) and the UV background will reduce the effectiveness of 
cooling in low density highly ionized gas (Efstathiou 1992). 
A UV background can therefore suppress the collapse of gas 
in regions of low overdensity. It is this low density photoion- 
ized gas that we believe accounts for the Lya absorption 
lines (Cen et al. 1994, Hernquist et al. 1996, Theuns et 
al. 1998, Bryan et al. 1999). Photoionization can also sup- 
press the collapse of gas in haloes with circular speeds of up 
to Vc ~ 20-30kms~^. However, numerical simulations have 
shown that a UV background cannot prevent the collapse of 
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gas in haloes with higher circular speeds, though it can re- 
duce significantly the efficiency with which low density gas is 
accreted onto massive galaxies (Quinn, Katz and Efstathiou 
1996, Navarro and Steinmetz 1997). 

To explain the galaxy luminosity function, feedback is 
required in galaxies with circular speeds Vc SL 50 kms^^ with 
characteristic virial temperatures of 10^ K. Energy injec- 
tion from supernovae is probably the most plausible feed- 
back mechanism for systems with such high virial tempera- 
tures. Winds from quasars might also disrupt galaxy forma- 
tion (Silk and Rccs 1998) or, more plausibly, limit the growth 
of the central black hole (Fabian 1999). Here we will be con- 
cerned exclusively with supernova driven feedback and will 
not consider feedback from an active nucleus. Simple para- 
metric models of supernovae feedback were developed by 
White and Rces (1978) and White and Frcnk (1991) and 
form a key ingredient of semi-analytic models of galaxy for- 
mation {e.g. Kauffmann, White and Guiderdoni 1993, Lacey 
et al. 1993, Cole et al. 1994, Baugh et al. 1996, 1998, 
Somerville and Primack 1999). In this paper, we develop a 
more detailed model of the feedback process itself. Previ- 
ous papers on supernovae feedback include those of Larson 
(1974), Dekel and Silk (1986) and Babul and Rees (1992). 
These authors compute the energy injected by supernovae 
into a uniform interstellar medium (ISM) and apply a sim- 
ple binding energy criterion to assess whether the ISM will 
be driven out of the galaxy. The feedback process in these 
models is explosive, operating on the characteristic timescale 
of ~ lO" lO'^ yrs for supernova remnants to overlap. This 
is much shorter than the typical infall timescale of hot gas 
in the halo, begging the question of how a resevoir of cold 
gas accumulated in the first place. The present paper differs 
in that we model the ISM as a two-phase medium consist- 
ing of cold clouds and a hot pressure confining medium, i.e. 
as a simplified version of the three-phase model of the ISM 
developed by McKee and Ostriker (1977, hereafter M077). 
The cold component contains most of the gaseous mass of 
the disc and is converted into a hot phase by thermal evap- 
oration in expanding supernovae renmaiits. In this type of 
model, the cold phase can be lost gradually in a galactic 
wind as it is slowly converted into a hot phase. 

The main result of this paper is that low rates of star 

formation can expel a large fraction of the baryonic mass 
in dwarf galaxies over a relatively long timescale of ~ 1 
Gyr. We therefore propose that effective feedback can oper- 
ate in an steady, unspectacular mode; strong bursts of star 
formation and superwind-like phenomena {e.g. Heckman, 
Armus and Miley 1990) are not required, although galax- 
ies may experience additional feedback of this sort. In fact, 
hydrodynamic simulations suggest that nuclear starbursts 
are ineffective in removing the ISM from galaxies with gas 
meisses ^ W^Mq (Mac Low and Ferrara 1999, Strickland 
and Stevens 1999) because hot gas generated in the nuclear 
regions is a expelled in a bipolar outflow without coupling 
to the cool gas in the rest of the disc. This result provides 
additional motivation for investigating a "quiescent" mode 
of feedback. Silk (1997) describes a model which is similar, 
in some respects, to the model described here. However, the 
model described here is more detailed and allows a crude 
investigation of the radial properties of a disc galaxy during 
formation. 



The layout of this paper is as follows. A simple model 
of star formation regulated by disc instabilities is described 
in Section 2. This is applied to 'closed box' {i.e. no infall or 
outflow of gas) models of disc galaxies neglecting feedback. 
Section 3 describes a model of the interaction of expanding 
supernovae shells in a two-phase ISM. This section is based 
on the model of M077, but instead of focussing on equi- 
librium solutions that might apply to our own Galaxy, we 
compute the net rate of conversion of cold gas to hot gas 
incorporating the model for self-regulating star forniatioii. 
This yields the temperature and density of the hot phase 
as a function of time and radius within the disc. Section 4 
revisits the model of Section 2, but includes simultaneous 
infall and outflow of gas. This model is extended in Section 
5 to include a galactic fountain, the pressure response of the 
cold ISM to the hot phase, and a model of chemical evolu- 
tion. Section 6 describes some results from this model and 
discusses the effects of varying some of the input param- 
eters. In addition, the efficiency of feedback is computed 
as a function of the circular speed of the surrounding dark 
matter halo. Our conclusions are summarized in Section 7. 
Although we focus on disc galaxies in this paper, a simi- 
lar formalism could be applied to the formation of bulges 
if the assumption that gas conserves its angular momentum 
during collapse is relaxed. 



2 STAR FORMATION REGULATED BY DISC 
INSTABILITIES 

2.1 Rotation curve for the Disc and Halo 

The dark halo is assumed to be described by the Navarro, 
Prenk and White (1996, hereafter NEW) profile 



pH{r) 



ScPc 



{cx){l + cx)^ 



''■/r-u, 



(2) 



where pc is the critical density, rv is the virial radius at 
which the halo has a mean overdensity of 200 with respect 
to the background and c is a concentration parameter (ap- 
proximately 10 for CDM models). The circular speed corre- 
sponding to this profile is 



(3) 



2 , s 2 1 [ln(l + ex) - cx/{l + cx)] 2 GM, 

VH{r) =v^- ' ; \ \ v^ = 

X [ln(l + cj - c/(l -I- c)] rv 

where M„ is the mass of the halo within the virial radius. 

We assume that the disc surface mass density distribu- 
tion is described by an exponential, 

P-D{r) = Mo exp(-r-/rD), Md = 2nfiorj), (4) 

where Md is the total disc mass. The rotation curve of a 
cold exponential disc is given by (Freeman 1970) 

v%{r) = 2vly' [Io{y)Ko{y) - h{v)Kr{y)] , 



(5) 



1 r 

2 m 



2 



GMe 



To relate the disc scale length, ro, to the virial radius 
of the halo rv, we assume that the angular momentum of 
the disc material acquired by tidal torques is approximately 
conserved during the collapse of the disc (see Fall and Efs- 
tathiou 1980). This fixes the collapse factor 



J coll — 

ro 



(6) 
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in terms of the dimensionless spin parameter Xh = 
J\E\^^^G~^ M~^^^ of the halo component. The spin param- 
eter is found to have a median value of ~ 0.05 from N-body 
simulations (Barnes and Efstathiou 1987), and for the mod- 
els described here, this value is reproduced for collapse fac- 
tors of around 50. A more detailed calculation of the collapse 
factor of the disc is given in Section 4. 

2.2 Vertical scale height of the disc 

The velocity dispersion of the cold gas clouds in the vertical 
direction is assumed to be constant and equal to ag. The 
equations of stellar hydrodynamics then give the following 
solution 



(7) 



where /ig is the surface mass density of the gas and the scale 
height is given by 



(8) 



Equation must be modified to take into account the steL 
lar disc. We do this approximately by assuming 'disc pres- 
sure equilibrium' (Talbot and Arnett 1975) 

„2 



H. 



1 



(9) 



" nGfig{l+l3/a)' 

where the quantities a and P relate the vertical velocity 
dispersion and surface mass density /i* of the stars to 
those of the gas clouds 



U ^ LtC g , 



(10a) 
(10b) 



2.3 Stability of a two-component rotating disc 

The stability of rotating discs of gas and coUisionless par- 
ticles to axisymmetric modes has been analysed in classic 
papers by Goldreich and Lynden-Bell (1965) and by Toomre 
(1964). Here we use the results of Jog and Solomon (1984) 
who analysed the stability of a rotating disc consisting of 
two isothermal fluids of sound speeds ci , C2 and surface mass 
densities /ii and fi2- These authors find that such a disc is 
stable to axisymmetric modes of wavenumber k if 



k 



+ 



2-kGh2 



< 1,(11) 



(12) 



Equation ( |ll| ) yields a cubic equation for the most unstable 
mode km ■ Solving this equation in terms of the parameters a 
and p of equations (10), and ignoring the small differences 
between a gaseous and coUisionless disc, we can write the 
stability criterion for a two-component system as 



K? (l-Ffc2cf/«:2) k2 (l+fc2c2/fi:2) 

where k is the epicyclic frequency 
K = 2a; l-| 



'-9{oc,P). 



(13) 



This is identical to the Goldreich-Lynden-Bell criterion ex- 
cept for the factor g{a, 0). This factor is plotted in Figure 1 
for various values of ct and p. 



a" 




B 10 



Figure 1. The factor g{a, P) appearing in the stability criterion 
of equation (13) plotted against /3 for three values of a 



2.4 Star formation and supernovae energy input 

We assume a stellar initial mass function (IMF) of the stan- 
dard Salpeter (1955) form 



dN, 
dm 



Am' 



mi < m < mu, x = 1.35, 



(14) 



mi = O.IMq, mn = 5OM0, 



and that each star of mass greater than 8Mq releases 
lO^^-Esi ergs in kinetic energy in a supernova explosion. For 
the IMF of equation (|l^, one supernova is formed for ev- 
ery 125M0 of star formation. The energy injection rate is 
therefore related to the star formation rate by 

Esn = 2.5 X W^^EsiAL erg/sec, (15) 

where M, is the star formation rate in Mq per year. 



2.5 Energy dissipated by cloud collisions 

We assume cold clouds of constant density = 7 x 10"^'^ 
g/cm"' with a distribution of cloud radii 

dNca 



da 



= Noa~ 



ai < a < a„, 



(16) 



ai = 0.5 pc, a„ = 10 pc, 



(M077). Following M077, the clouds are assumed to have 
an isotropic Gaussian velocity distribution with velocity dis- 
persion independent of cloud size and that they lose energy 
through inelastic collisions. The rate of energy loss per unit 
volume is given by 



dE, 



coll 



dtdV '^'i^^^'^PcNcldfo-gla, 
a^/ai facial (^^^)2 

(x^ + j/-^) X y ' 



(17) 



where Nd is the local cloud density Nd = No /3a; . Inte- 
grating equation (17) over the vertical direction and using 
equation (W for the scale height, the rate of energy loss per 
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Table 1: Parameters of Model Galaxies 



Vc (km/s) JJmax (km/s) Vy/Vc 

Model MW 280 212 0.45 

Model DW 70 53 0.45 



ro (kpc) Md (Mq) fcoii c Xh 
3.0 5.5 X lO^o 50 10 0.065 

0.2 2.3 X 10* 50 10 0.065 



unit surface area E^ou is 

eZu = 5.0 X 10"'' (l + ^) fi%a,5 erg/sec/pc^ (18) 

where figs is the surface mass density of the gas component 
in units of 5Mq/pc^ and CTgs is the cloud velocity dispersion 
in units of 5 km/sec. These values are close to those observed 
in the local solar neighbourhood. To estimate the efficiency 
with which supernovae accelerate the system of clouds, we 
normalize to the observed net star formation rate of the 
Milky Way. Assuming that the gas distribution has a flat 
surface mass density profile to Rmax ~ 14 kpc (Mihalas and 
Binney 1981), /3 « 10, a ~ 5, and equating the integral of 
( p^ to tciSsn (equation p^) , we find 

ecEsiM, = 0.004. (19) 

An efficiency parameter of ec = 0.01 produces a net star 
formation rate of O.4M0 /yr which is reasonable for a Milky 
Way-like galaxy. We will therefore adopt a constant value of 
Ec ~ 0.01 in the models of the next subsection. The value 
of Ec will, of course, depend on the properties of the clouds, 
ISM and star formation rate. For example, in the model of 
M077 the clouds are accelerated by interactions with the 
cold shells surrounding supernova remnants and they find 
efficiencies Ec of typically a few percent. We investigate the 
effect of varying Ec in Section 6. 



2.6 Self-regulating models without inflow or 
outflow 

The equations derived above allow us to evolve an initially 
gaseous disc and to compute the local star formation rate, 
cloud velocity dispersion etc. The system of stars and gas is 
constrained to satisfy the stability criterion of equation (p^, 
which fixes the cloud velocity dispersion ag. There is some 
empirical evidence that star formation in nearby galaxies 
is regulated by a stability criterion of this sort {e.g. Ken- 
nicutt 1998). The energy lost in cloud collisions (equation 
p^ ) is balanced against the energy input from supernovae 
assuming a constant efficiency factor Ec = 0.01. We assume 
further that a — 5 (equation 10a) i.e. that stars are in- 
stantaneously accelerated to higher random velocities than 
the system of gas clouds, and that the properties of the gas 
clouds (mass spectrum, internal density, etc) are indepen- 
dent of time. These are clearly restrictive assumptions, but 
they allow us to generate simple models of self-regulating 
star formation with only one free parameter Ec. 

We study the evolution of two model galaxies with 
parameters listed in Table 1. Model MW has parameters 
roughly similar to those of the Milky Way and model DW 
has parameters similar to those of a relatively high surface 
brightness dwarf galaxy. 



Figure 2 shows the evolution of the gas and stellar sur- 
face mass densities of the two models. The net star forma- 
tion rates, gas fractions and mean gas cloud velocity disper- 
sion are plotted in Figure 3. In model MW, the star forma- 
tion rates are initially high (> lOOM0/yr) and hence the 
timescale for star formation is short; half the disc mass is 
converted into stars in 10^ years. The star formation rate de- 
clines rapidly to less than lM0/yr after a few Gyr. As figure 
2 shows, the star formation at early times is concentrated 
to the inner parts of the disc which have a high surface den- 
sity and hence the gas distribution develops a characteristic 
surface density profile with an inner 'hole', similar to what 
is seen in the HI distributions in real galaxies (see Burton 
1991). The stellar disc is truncated at about the Holmberg 
(1958) radius (r/ro « 5), in rough agreement with obser- 
vations. The truncation arises because the gas disc becomes 
thick at large radii (equation ^) and the rate of energy lost in 
cloud collisions can be balanced by a very low star formation 
rate. 

The evolution of model DW is qualitatively similar, 
though the star formation rate is scaled down roughly in 
proportion to the disc mass. Half the gas is converted to 
stars by 3 X 10^ yr, and the gas fraction is 0.12 after 10^° 
yr, similar to the final gas fraction of 0.13 in model MW. 

Neither of these models is satisfactory. The star for- 
mation rate in model MW is too high at early times to be 
compatible with deep number counts (see e.g. Ellis 1997), 
which require more gentle star formation rates in typical L* 
galaxies. Model SW converts most of its gas into stars on 
a short timescale and so does not solve the problem raised 
in the introduction of explaining the flat faint end slope of 
the luminosity function in CDM-like models. As we will see 
in later sections, infall of gas provides the solution to the 
former problem, since this allows the disc to build up grad- 
ually on a cooling or dynamical timescale. Outflow of hot 
gas heated by supernovae provides a solution to the latter 
problem. 



3 EVOLUTION OF A TWO-PHASE ISM 

In this Section we consider the interaction of a multiphase 
interstellar medium with expanding supernova remnants fol- 
lowing the model of M077 and discuss the conditions under 
which a protogalaxy can form a wind. The key ingredients 
of the model are as follows. Most of the cold gaseous mass 
is assumed to be in cold clouds with properties as given in 
Section 2.5. Supernovae explode and their remnants propa- 
gate evaporating some of the cold clouds and forming a low 
density hot phase of the ISM. The star formation rate there- 
fore determines the evaporation rate and hence the rate of 
conversion of the cold phase to a hot phase. A wind from the 
galaxy can result if the hot phase is: (i) sufficiently pervasive 
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Figure 2. The evolution of the gas (solid lines) and stellar (dashed lines) surface mass density distributions according to the simple 
self-regulating model described in this section. The results are shown for ages of 0, 0.1, 1, 3, 6 and 10 Gyr. 
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Figure 3. Evolution of the star formation rate, gas fraction and gas cloud velocity dispersion in the self-regulating model. 
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(filling factor of order unity) , (ii) low density (so that radia- 
tive cooling is unimportant) and (iii) the temperature of the 
hot phase exceeds the virial temperature of the galaxy. In 
this section we follow closely the theory of the ISM devel- 
oped by M077 and we use their notation where possible. 



3.1 Evaporation of cold clouds 

An expanding supernovae remnant will evaporate a mass of 



Me„ « 540 Egi^S 



6/5^,-3/5^-4/5 



Me, 



(20) 



where Uh is the density interior to the supernovae remnant 
and E (in pc^) is the evaporation parameter introduced by 
M077 



E = 



7 



AnaiNcirtiK 



(21) 



Here the parameter 7 relates the blast wave velocity to the 
isothermal sound speed (nj, — ych, 7 ~ 2.5) and the pa- 
rameter c/>re quantifies the effectiveness of the classical ther- 
mal conductivity of the clouds (kb// = k0k) and so is less 
than unity if the conductivity is reduced by tangled mag- 
netic fields, turbulence etc. Using equations (j^ and (^ to 
estimate the mean cloud density we find 



280-2 



95 



1 



M?5 (1 + 



pc = /eEc 



95 pc^(22) 



where E© is the evaporation parameter characteristic of the 
local solar neighbourhood (P/a ~ 2). 
Evaluating equation (EOl), we find 



Me„ ^ 1390 Eli^f-^'^^Tnlt'^ Mq, 



(23) 



where nh-2 is in units of 10~^cm~^ (a characteristic value 
for the hot component). Thus, provided thermal conduction 
is not highly suppressed, a single supernovae remnant can 
evaporate a much larger mass than the 125 AIq formed in 
stars per supernovae for a standard Saltpeter IMF (Section 
2.4). If a significant fraction of this evaporated gas can es- 
cape in a wind, then star formation will be efficiently sup- 
pressed. 



3.2 Temperature and density of the hot phase 

To compute the properties of the hot phase we assume that 
the disc achieves a state in which the porosity parameter Q 
is equal to unity. The disc is then permeated by a network of 
overlapping supernovae remnants. Ignoring cooling interior 
to the remnants (which we will see is a reasonable approxi- 
mation for an ISM with low metallicity) the age, radius and 
temperature of a SNR when Q = 1 are given by 



to = 5.5 X lo'^s; 

Ro = lOOS- 



5/11^,-6/11 p-3/ll 3/11 



7 



13 



2/11^,2/11^1/11^-1/11 



7 



pc. 



T„ = 1.2 X 10*5,«/"7-«/"i5!/"n:«/" K. 



(24a) 
(24b) 
(24c) 



^51 '-h 

where S13 is the supernova rate in units of 10~'^^pc~^yr~^. 
The density of a remnant at to (n^ ~ Met,/(4/37rJ?o)), gives 
an approximate estimate of the density of the ambient hot 
phase 



rih ^ 4.3 X 10 5i3 7 /j. 



(25) 



Inserting this estimate into equations (24) we find 

J- 1 o„ in6c<-0-36 -0. 64 /17, r -,-0.11 

to = 1.2 X 10 7 (-E51./E) yr. 



(26a) 
(26b) 

To = 6.6 X 10"(Si3£;5i/e/7)°-'' K, (26c) 



D laA/Q I \ -0.21 17,0.04 _fO. 035 

Ro = 164(513/7) £51 /e pc 



and the rate at which clouds are evaporated is 

A> r, -7 w 1 n-10 cO.71 0.29 ^0.71 r-0. 29 -3 -1 /o-tN 

Me„ = 2.7 X 10 Sn 7 h Mopc yr . (27) 

Integrating equation (^^ over the scale height of the disc 
gives the evaporated mass per unit area, 



1 X 10"^ 



'95 



Ms5 (!+/?/«) 



crO.71 0.29 EiO.Tl /.-0.29 T.r -2 -1 

-513 7 -E51 h M0PC yr . 



(28) 



Adopting a cooling rate of A 2.5 x 10 Tj erg cm^ 



for 10" 



. r < 10"= for a gas with primordial composition. 



the ratio of to to the cooling time tcooi is 

to 



tc 



o.5r5 /e , 



(29) 



Thus if the temperature of the hot phase is higher than 
about 10"K, the assumption that cooling can be neglected 
will be valid. A cooling function for a gas with primor- 
dial composition will be used throughout this paper. As the 
metallicity of the gas builds up, the cooling time of the hot 
component will shorten and more of the supernovae energy 
will be lost radiatively. This effect will reduce the efficiency 
of feedback in galaxies with high metallicity but is not in- 
cluded in this paper. 



3.3 Simple self-regulating model with outflow 

In this section we apply the results of the previous para- 
graphs to construct a simplified self-regulating model with 
outflow. The star formation rate is governed by the self- 
regulation algorithm as in Section 2.6 with the parameter 
Ec = 0.01. This provides an estimate of the local supernova 
rate per unit volume which we insert in equations (26) to 
compute the properties of the hot phase, adopting a value 
<;/>« = 0.1 in equation ( pl|) for the conduction efficiency pa- 
rameter. The hot gas will be lost from the system if its spe- 
cific enthalpy 

1 2 5p 
-V H 

2 2p 

exceeds to within a factor of order unity its gravitational 
binding energy per unit mass. If the gas has an initial 
isothermal sound speed of d — (kT/np) = 37T5^'^kms~^ 
(for a mean mass per particle of fip = 0.61mp), conserva- 
tion of specific enthalpy implies that the wind will reach a 
bulk speed of « V^d. Some of the thermal energy will 
be lost radiatively, and in fact the spherical steady wind 
solutions described in Appendix B suggest that a more ac- 
curate criterion for the wind to escape from a galaxy is 
Vuj ~ \/2.5ci > Vesc, where Vesc is the escape speed from 
the centre of the halo (neglecting the potential of the disc). 
If \/2.5ci < Vesc, the hot phase is returned instantaneously 
to the cold phase. This type of binding energy criterion for 
outflow has been adopted in previous studies {e.g. Larson 
1974, Dekel and Silk 1986) and is clearly oversimplified, as 
are the assumptions of instantaneous mass loss and return 
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Figure 4. The left hand panels show the evolution of the gas (solid lines) and stellar (dashed lines) surface mass density distributions for 
ages of 0, 0.1, 1, 3, 6 and 10 Gyr as in figure 2. The panels to the right show various properties of the hot gas component as a function of 
the disc radius r/r£). The solid lines show the temperature, dashed lines show the density and the dotted lines show the ratio of overlap 
to cooling timescales, to/tcool- 



of cold gas. These points will be discussed further in Section 
6, but for the moment these assumptions will be adopted 
to illustrate the qualitative features of the model. As gas is 
lost from the system, the circular speed of the disc compo- 
nent (equation 5) is simply rescaled by the square root of 
the mass of the disc that remains. 

The evolution of the surface mass densities for the two 
disc models is illustrated in Figures 4 and 5. In model MW, 
the evolution is similar to that without outflow shown in 
Figure 0. With the simple prescription for mass loss used 



here, no hot gas is lost unless the temperature of the hot 
phase exceeds Tcrit « 5 x lO^K. This does happen at early 
times when the star formation rate is high, and about 25% 
of the galaxy mass is lost within 10^ yr. Thereafter, no more 
mass is lost and a nearly exponential disc is built up with 
a gas distribution containing a central hole as in Figure 2. 
The star formation rate in this model declines strongly with 
time, exceeding lOOA/0/yr in the early phases of evolution. 

The behaviour of model DW is qualitatively different. 
Here the critical temperature for mass loss is much lower. 
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(a) Model MW 




t (Gyr) 

Figure 5. Evolution of the star formation rate, gas fraction and 

Tcrit ~ 3 X 10^ K, hence half the mass of the galaxy is 
expelled by ^ 10* yr and and 66% by 1 Gyr. After 1 Gyr, 
the temperature of the hot phase drops below Tcrit and the 
galaxy settles into a stable state with a low rate of star 
formation. 

The wind prescription in these models, and particularly 
the assumption that gas below the critical temperature nec- 
essary for escape is returned instantaneously to the cold 
phase, is clearly oversimplified and so the mass loss fractions 
should not be taken too seriously. A more detailed model 
is developed in Section 6. A more serious deficiency of the 
model presented here is that the entire gas disc is assumed 
to have formed instantaneously ai t — 0. This is unrealistic 
and leads to high rates of star formation and gas ejection at 
early times. A simple infall model, similar to those adopted 
in semi-analytic models (White and Frenk 1991, Cole et al. 
1994) is included in the next Section. 



4 INFALL MODEL 

4.1 Conservation of specific angular momentum 

Following Fall and Efstathiou (1980), the gas is assumed to 
follow the spatial distribution of the halo component with 
the same distribution of specific angular momentum prior 
to collapse. The halo is assumed to rotate cylindrically with 
rotation speed v^'^{zajh), where ujh is the radial coordinate 



(b) Model DW 



1 r 




i — I I I I ll [ 



Q L I I I I I I I I 

0.01 0.1 1 10 

t (Gyr) 

gas cloud velocity dispersion for the models shown in figure 4. 

in the cylindrical coordinate system. The gas is assumed to 
conserve its specific angular momentum during its collapse, 
so that the final specific angular momentum of the disc at 
radius zud, ho = vjdv^q^, is equal to the specific angular 
momentum of the halo ho ~ ■cohv^h^ at the radius vjh from 
which the gas originated. Mass conservation relates the radii 
zuh and zud, 

dujH _ M-p(^p) Mh -too ^gp^ 

where Mh /Md is the ratio of the halo to disc mass interior 
to the maximum infall radius of the disc (see figure 6a below) 
and is the projected surface mass density of the halo 

f oo 

AtH(tJ7)=2/ ph{{-oo'' + zy^)dz. (31) 
Jq 

The solution of equation ( po| ) yields ■cod{'coh) and the ro- 
tation speed of the halo follows from the conservation of 
specific angular momentum, wjf* = vj dv^d^ {vj d) / tjo h ■ The 
results for the parameters of models MW and DW are shown 
in figure where we have used the notation s — vj/rn- 
When expressed in the dimensionless units of Figure ^ the 
solutions for models MW and DW are identical. 

This prescription is guaranteed to form an exponential 
disc with the required parameters. The derived rotation ve- 
locity of the halo is almost independent of radius in general 
agreement with what is found in N-body simulations (Frenk 
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0.01 0.1 1 10 100 10 20 30 40 

Figure 6. Figure 6a shows the solution of tlic mass conservation equation (30) relating the initial halo radius sh = t^h/td to the 
final disc radius so = ^d/'>'d- The solid line in figure 6b shows the derived rotation curve of the halo in units of Vc = (GMjy/ro)^^^ 
assuming conservation of specific angular momentum hn = ho. The dashed line shows the fitting function of equation (39). 




Figure 7. Tiic free-fall (solid lines) and cooling times (dashed lines) for the two model galaxies plotted as a function of halo radius 
^h/^d- Note that for these models, the ratio of baryonic to dark mass within the virial radius is 0.1. 

et al. 1988, Warren et al. 1992). The upturn in the halo ro- that the the halo rotation velocity is constant at O.OQSwc at 

tation speed at snimax) « 30 is caused by the rapid decline large radii, 
in the mass of the input exponential disc at large radii and 
is of little consequence in the discussion that follows. The 
values of the spin parameter quoted in Table 1 were derived 
from the mass and binding energy of the halo and assuming 
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4.2 Mass infall rate 

To determine the gas infall rate we compute the free fall 
time for a gas element at rest at radius n, 



ts 



dr 



V2[<^i/(rO-0H(r)]i/2' 
and the cooling time 
_ 3 fcT„ X 1.92 

^cool — 7T 



(32) 



(33) 



2 A(T„)ne(r)' 

where n^jr) is the electron density. The temperature T„ in 



equation ( p3D is set to the virial temperature derived from 
the equation of hydrostatic equilibrium assuming that the 
temperature is slowly varying with radius 



Tv ~ -VH{r)- 



k d\npb{ry ^"^^^ 

where we assume that the baryons follow the same spatial 
distribution as the halo. The infall rate is given by 



Minf = 'inpb(rH)rH< 



dvH (iff = t)/dt 



iff > icooi 



drnitc 



t) jdt 



tc 



.(35) 



Finally, conservation of specific angular momentum specifies 
the final radius in the disc for each gas element. Since the 
halo is assumed to rotate on cylinders, the gas near to the 
poles in an infalling shell has a lower specific angular mo- 
mentum that the gas at the equator. The infalling material 
is therefore distributed through the disc according to 

2TT-UJ DflD{'!^D)duJD = ^inf- ^„2 ^ \l/2 ' ^'^^'^ 



where vjd and vjh are related by the solution of equation 

Equations (B2|) - ( |36[ ) specify the infall model. The free- 
fall and cooling times of the two model galaxies are shown 
in Figure 0. In the larger galaxy, gas within rn/rD ~ 10 
infalls on the free-fall timescale and ends up within one scale 
length of the final disc. The material in the outer parts of the 
disc infalls on the cooling timescale. In contrast, apart from 
a small amount of gas in the very central part of the halo 
with virial temperature < 10* K, the gas in the dwarf galaxy 
infalls on a free-fall timescale because the cooling time is so 
short. 



4.3 Simple self-regulating model with inflow and 
outflow 

The models described in this section are exactly the same 
as those described in section 3.3, except that we grow the 
discs gradually using the infall model of sections 4.1 and 
4.2. In the models described below, inflow and outflow are 
assumed to occur simultaneously. This is often assumed in 
semi-analytic models of galaxy formation {e.g. Cole et al. 
1994, Somerville and Primack 1999) and may not be com- 
pletely unrealistic if the infalling gas is clumpy. The dark 
matter haloes will contain significant sub-structure {e.g. 
Moore et al. 1999) which may contain pockets of cooled 
gas. Furthermore, if the cooling time is short compared to 
the dynamical time, the infalling gas will be thermally un- 
stable (Fall and Rees 1985) and will fragment into clouds. 
These will fall to the centre on a free-fall timescale if they 



are sufficiently dense and massive that gravity dominates 
over the ram pressure of the wind. This requires clouds with 

masses 



mcioud Z 9.5 X 10^ M( 



(^clouc 
Ikpc 



Ikpc / V lOkpc 



VlMo/yr/ UOOkm/s/ UOOkm/s/ 



(37) 



where flcioud is the radius of the cloud. However, even if 
( p7| ) is satisfied, the clouds may be sheared and disrupted 
into smaller clouds by Kelvin-Helmholtz instabilities on a 
timescale of a few sound crossing times as they fiow through 
the wind {e.g. Murray et al. 1993). The wind energy will 
be partially thermalized in shocks with the infalling clouds 
and dissipated in evaporating small clouds. But for the 
typical mass outflow rates expected from dwarf galaxies 
{Mw ^ 0.2Mq /yi), the rate at which energy is supplied by 
the wind Em = l/2Mn,v^ is much smaller than the energy 
lost in radiative cooling. 



—. — ^ 50A_23 



yiOOkm/s/ VlOkpc 



Mm 



yiOOkm/sJ yiAfo/yr 



(38) 



where rcooi is the radius at which the cooling time is equal 
to the age of the system. 

The qualitative picture that we propose is as follows. In 
galaxies with a short cooling time, clouds formed by thermal 
instabilities will infall ballistically if (^^ is satisfied. If (js^ 
is not satisfied, the ram pressure of the wind will drive out 
the infalling gas and infall will be suppressed. With infall 
suppressed, the star formation rate in the disc and the wind 
energy will decline until infall can begin again. The wind will 
be partially thermalised before reaching rcooi and completely 
thermalised at ~ Tcooi, but the energy supplied by the wind 
will be small compared to the energy radiated by the gas 
at r ^ '"cool and so cannot prevent radiative cooling. If ( |37[ ) 
is satisfied, some of the outfiowing gas may fall back down 
to the disc after shocking against infalling clouds. However, 
in the models described here the efficiency of converting in- 
falling gas into stars is low in dwarf systems, so provided the 
gas does not cycle around the halo many times, neglecting 
return of some of the outfiowing gas should not afi'ect the 
qualitative features of the models. The global geometry of 
the system, e.g. if the wind is weakly coUimated perpendic- 
ular to the disc, may also permit simultaneous infiow and 
outflow of gas. 

The interaction of an outflowing wind with an inhomo- 
geneous infalling gas clearly poses a complex physical prob- 
lem. In reality, the process may be far from steady, with out- 
flow occurring in bursts accompanied by infall from discrete 
sub-clumps containing cooled gas. In the models described 
below and in the rest of this paper, we will assume that the 
infall and outflow occur simultaneously, steadily and with- 
out any interaction between the inflowing and outflowing 
gas. As the discussion of the preceeding two paragraphs in- 
dicates, this is obviously an over-simplification. It should be 
viewed as an idealization, on a similar footing to some of 
the other assumptions adopted in this paper {e.g. spherical 
symmetry, neglect of halo substructure and merging, steady 
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Figure 8. The left hand panels show the evolution of the gas (solid lines) and stellar (dashed lines) surface mass density distributions 
for ages of 0, 0.1, 1, 3, 6 and 10 Gyr as in figure 2. The panels to the right show the radial distribution of density, temperature and ratio 
of overlap to cooling timescales for the hot gas component. 



star formation rates etc) designed to give some insight into 
how a quiescent mode of feedback might operate. 

The analogues of figures ^ and |^ for the models incor- 
porating infall and outflow are shown in figures ^ and ^. The 
discs build up from the inside out, as in the models of disc 
formation described by Fall and Efstathiou (1980) and Gunn 
(1982). Most of the star formation occurs in a propagating 
ring containing the most recently accreted gas. The most sig- 
nificant differences from the models of section 3.3, are the 
net rates of star formation (figure ^) and the timescale of 
outflow. The initial high rates of star formation in the mod- 
els of section 3.3 are suppressed in the models with infall, 
and the timescale for outflow is now much longer because it 



is closely linked to the gas infall timescale. Apart from these 
differences, the flnal states, gas fractions and mass-loss frac- 
tions are similar to those in the models without infall. In 
model MW some outflow occurs when t lO^yrs and the 
temperature of the hot gas is high enough that it can escape 
from the system. Thereafter, the hot component cannot es- 
cape and the disc builds up without further outflow. About 
17% of the total galaxy mass is expelled in the early phases 
of evolution, but as we have described above, this could be 
an overestimate since some of this gas may be returned to 
the galaxy if the wind energy is thermalized before it reaches 
the virial radius. In contrast, model DW drives a wind until 
t ^ 1 Gyr and expels about 74% of its mass. About half of 
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(a) Model MW (b) Model DW 




the gas is lost within 3 x 10 yrs, i.e. on about the infall 
timescale for most of the gas in the halo (see figure 7b). 



5 REFINEMENTS OF THE MODEL 

The models described in the previous sections contain a 
number of simplifications, which we will attempt to refine 
in this section. We do not address the problem of the inter- 
action of a wind with the infalling gas, which is well beyond 
the scope of this paper. Instead, we introduce some simple 
improvements to the infall model (§5.1), model for mass loss 
from the galactic disc (§5.2, §5.3) and the pressure response 
of the cold ISM to the hot phase (§5.4, §5.5). 

5.1 Infall Model 

In Section 4, we used a simplified model of infall that guar- 
antees the formation of an exponential disc if angular mo- 
mentum is conserved during collapse. In this section, we 
assume a specific functional form for the rotation velocity of 
the halo, 

rot/ \ \^/^'^} / /on\ 

vh [s) = ciVc ,^ , , , , , , . , s = rrD (39) 

(1 + (s/C2))(l -I- (s/C3)=4) 

with ci — 0.115, C2 = 0.6, C3 = 16, C4 = 0.25. The functional 
form and coefficients in equation ( ^9| ) have been chosen to 
provide a good fit to the halo rotation profile derived in 
section 4.1 from conservation of specific angular momentum 



and is plotted as the dashed line in figure (6b). As in the 
previous section, the gas is assumed to follow the same ra- 
dial density distribution and rotation velocity as the halo 
component, but its final radius in the disc is computed by 
assuming conservation of specific angular momentum and 
self-consistently solving for the rotation speed of the disc 
component. The halo component is assumed to be rigid and 
the contribution of the disc component to ujf' is computed 
using the Fourier-Bessel theorem (see Binney and Tremaine 
1987 §2.6) 

v^{r) = -r / S{k)Ji{kr)k dk, 
Jo 

/> oc 

S{k)^-2nG Jo{kr)fiD{r)r dr. (40) 
^0 

Equation ( |4C| ) is time consuming to evaluate accurately and 
in our application «^(r) must be computed many times. A 
fast algorithm has therefore been developed as described in 
Appendix A. The epicyclic frequency k is required in equa- 
tion (|l^ to compute the instantaneous star formation rate 
and is derived by numerically differentiating the rotation 
speed. 

With this formulation of the infall model, the in- 
fall rate is governed by the dark matter profile and the 
ratio of dark to baryonic mass within the virial radius, 
Mv/Md = (vv/vc)^ fcoii- For the models described here, we 
adopt M^/Md ~ 10, consistent with the parameters listed 
in Table 1. The final disc surface mass density will be close 
to an exponential by construction, since the halo rotation 
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velocity (|39|) has been chosen to match the rotation profile 
derived by assuming an exponential disc and conservation 
of specific angular momentum. 



5.2 Galactic Fountain 

In previous sections, we have assumed that gas is lost from 
the disc if the bulk velocity of the wind Hn, ~ \/2.5ci, ex- 
ceeds the escape speed Vesc from the halo, but is otherwise 
returned instantaneously to the ISM. More realistically, gas 
with < Vesc will circulate in the halo along a roughly 
ballistic trajectory and will cool forming a galactic fountain 
(Shapiro and Field 1976, Bregman 1980). In the models de- 
scribed in this section, hot gas with Vm < v^sc is returned 
to the disc at the radius from which it was expelled after a 
time tret, 



(41) 



i.e. we ignore the gravity of the disc and the angular mo- 
mentum of the gas in computing the ballistic trajectory of 
a gas element. 



5.3 Escape Velocity of the Wind 

The detailed dynamics of the hot corona itself is complicated 
and beyond the scope of this paper. Type II supernovae at 
the upper and lower edges of the gas disc will be able to 
inject their energy directly into the hot gas, as will Type la 
supernovae forming in the thicker stellar disc. In addition, 
the hot component will interact with the primordial infalling 
gas in a complicated way as sketched in §4.3. 

In the absence of radiative cooling, the hot gas will ex- 
tend high above the galactic disc in an extended corona. For 
an isothermal corona, the equation of hydrostatic equilib- 
rium in the z-direction has the following approximate solu- 
tion, 



pjf (ro, z) 
Pff(O) 



: sech 



2P9 



— I sech 



2ps 



ii) 



exp 



zvu{r) 



dz 



2 , 2 
■OJ + Z , 



(42) 



where we have assumed that both the stars and the gas 
follow sech'^ vertical distributions (equation ^ and Ci is the 
isothermal sound speed of the hot gas. We define a character- 
istic scale height for the hot component, Hhoti''^), at which 
the density drops by a factor ~ e according to equation (^^. 
If radiative cooling were negligible, we would expect a sonic 
point in the flow at about z ~ Hhot. It is interesting to 
compare some characteristic numbers for the coronal gas: 



E^r^ ^ 1.1 X lO^^TfteMe. crg S-' 
-BsNii ^ 5.7 X 10'"'esNii£5iA^* erg s"' 

Hhot \ ( Rhot 



Eeool « 2 X 10^^nL2A-23 



1 kpc / V 3 kpc 



erg s 



(43a) 
(43b) 

"'(43c) 



where the rates Mev, Mt are in Mq /yr. Here Ein is the 
thermal energy injected into the hot corona by evaporating 



cold gas at a rate Mev, Esmi is the energy supplied to the 
corona by Type II supernovae forming above and below one 
scale height of the cold gas layer and the parameter esNii 
expresses the efficiency with which this energy is coupled 
to the hot coronal gas. -Bcooi is the rate of energy lost by 
a uniform density isothermal corona of scale height H^ot 
within a cylinder of radius Rhot. For a large galaxy such 
as the Milky Way that can sustain an evaporation rate of 
~ lOM0/yr, -Ecooi is small compared to Ei„ and it is a 
good approximation to neglect radiative cooling in the early 
stages of the flow (see Appendix B). However, in a dwarf 
galaxy -Ecooi is typically larger than Ein. In this case, we 
expect that the hot component will develop a sonic point at 
a characteristic cooling scale height -ffcooi 



-f^cool ^ ^uj^cool "^11 



100 km/s 



ThsK^.AZla kpc. (44) 



(see e.g. Kahn 1981 and Appendix B) and that most of the 
thermal energy will be converted into kinetic energy by the 
time the gas flows to i/cooi. We therefore ignore radiative 
cooling and estimate the bulk velocity of the wind at each 
radial shell in the disc from 



■^-^oLt^w = E2i + -EsNII, 

and to close the equations we assume that 

Mnii.t = = Me^>. 



(45a) 



(45b) 



(Note that the numerical coefficient in equation (43a) has 
been adjusted to give = \/2.5ci if -Esnii = so that the 
criterion for the wind to escape, > Vesc, is the same as 
in the preceeding Section). 

Energy input from Type II supernovae exploding high 
above the cold gas layer will make a small contribution to 
the thermal energy of the hot coronal gas. For values of 
E'J Pfjj will be about 20% or so of Ein and can- 
not be higher because Mf^ and the star formation rate are 
nearly proportional to each other (equation . Type la su- 
pernovae will also supply energy to the corona, with a time 
lag of perhaps ^ 1 Gyr (Madau, Delia Valle and Panagia 
1998). However, this effect will also make a small perturba- 
tion to the energy budget of the corona and so it is ignored 
here. Furthermore, in the models described here much of the 
gas is expelled on a timescale of ^ 1 Gyr, thus feedback is 
likely to be more or less complete before energy injection 
from Type la supernovae becomes significant. 



5.4 Pressure equilibrium and cold cloud radii 

In the models of Section 3 and 4, the cold cloud radii were 
kept constant irrespective of the pressure of the confining 
hot phase. More realistically, the cold cloud radii will adjust 
to maintain approximate pressure equilibrium with the hot 
phase. Thus 

where Uq denotes the cloud radii a,t the solar neighbourhood 



with values as given in equation ( |17[ ) and Tc is the internal 
temperature of the cold clouds. In our own Galaxy, pho- 
toelectric heating of dust grains is believed to be the main 
heating mechanism of the cold clouds (see e.g. Wolfire et al. 
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1995) but other heating mechanisms may be important, for 
example, cosmic-ray heating (Field, Goldsmith and Habing 
1969). We therefore expect that Tc varies in a complex (and 
uncomputable) way as a galaxy evolves. To assess the ef- 
fects of the pressure response of the cold clouds, Tc will be 
assumed to remain constant at 80K. The cloud radii are then 
determined solely by the pressure of the hot component via 
equation (|46[). 

The energy lost through cloud collisions (equation 
varies as (for fixed cloud masses). However, the cloud 
heating efficiency factor ec will also change as the cloud radii 
change in response to the pressure of the hot phase. In the 
model of M077, the energy acquired from momentum ex- 
change with cooling supernovae shells varies as a*. The net 
effect of these variations in the self-regulating star formation 
model is to introduce positive feedback, since a higher rate 
of star formation is required to balance energy lost through 
cloud coUisions in regions where the ambient pressure is 
higher. This is modelled by assuming EcoU oc {a/aQ)^ and 
Ec ~ ecoia/ uq)*' , where tc© is a fiducial efficiency factor. 



5.5 Induced star formation 

The maximum mass for an isothermal cloud in pressure equi- 
librium with the confining medium of pressure ph is given 
by the Bonner-Ebert criterion (Bonner 1956, Ebert 1955, 
Spitzer 1968), 



galactic fountain of metallicity Zp accreted at a rate dup/dt. 
The equation of chemical evolution is then 



^-3/2p-l/2 



rriBE = 1.18 

= 433 (^)^n,_2r,6)-'/'M0. 



(47) 



For our own Galaxy (M077), nn ~ 1.5 x 10-^cm-^ Tn ~ 
4 X 10^ and Ph ~ 3 x 10~^^dyne cm~^, hence ttibe ~ 
1700Mq. This is reasonably close to the upper mass limit, 
ruu = 43OOM0 , for the cold cloud mass spectrum adopted in 
this paper («„ = lOpc with pc = 7 x 10^^'^ g/cva?). Gravita- 
tional stability requires m„ ~ meE and we will henceforth 
impose this condition in determining the upper mass limit of 
the cold cloud spectrum. An increase in the pressure of the 
hot phase will lead to a decrease in m„ and hence to some 
pressure induced star formation. If the over-pressured clouds 
fragment into stars with an efficiency ebe, the induced star 
formation rate is given by 



dt 



Ml. 



2phln(m„/mL) | 



Oh/dt dph/dt > 
dph/dt < 



(48) 



where is the lower limit to the cloud mass spectrum, 
rriL 0.5Mq. This provides an additional source of positive 
feedback, since as the pressure of the hot component rises the 
star formation rate of the self-regulating model is enhanced 
by pressure induced star formation. 



5.6 Chemical evolution 

It is straightforward to include chemical evolution in the 
models using the instantaneous recycling approximation. We 
distinguish between 'primordial' infalling gas accreting at a 
rate djii/dt with metallicity Zi, and processed gas from the 



HgdZ = pd/is + [Zi — Z)djj.i + [Zf — Z)dfj,F, 



(49) 



(see e.g. Pagel 1997), where p is the yield. We adopt a yield 
of p = 0.02 and assume that the primordial gas has zero 
metallicity (Zj = 0). Gas ejected in a galactic fountain is 
assumed to have the same metallicity as the ISM at the 
time that the gas was ejected. Within the disc, the ISM gas 
is assumed to be perfectly mixed at all times. We normalize 
the metallicities to the solar value, for which we adopt Zq — 
0.02. 



6 RESULTS AND DISCUSSION 

6.1 Variation of input parameters 

In addition to the many simplifying assumptions introduced 
in previous sections, the model described here has 4 key 
parameters: (i) 0^, determining the efficiency of heat con- 
duction (equation GQ); (ii) tc©, controlling the star forma- 
tion rate (equation |19[); (iii) esNii, determining the efficiency 
with which energy fro m T ype II supernovae couples directly 
to the gas (equation 45a) ; (iv) ebe, setting the efficiency 



with which over- pressured ISM clouds collapse to make stars 
(equation ^). In addition, the ISM cloud radii can be al- 
lowed to vary in response to the pressure of the ISM as 
described in Section 5.4. 

The effects of varying these parameters are summarised 
in Table 2. Here we have run six models of galaxies MW and 
DW varying the input parameters. We list the final stellar 
mass Ms , gaseous disc mass Mg , and ejected mass M^j after 
10 Gyr for model DW (there is very little evolution after this 
time) and after 15 Gyr for model MW. The parameters fej 
and /* are the final ejected and stellar masses divided by the 
total baryonic mass (A^ej ~1" M^ -t- Mg^. T^j is the time when 
half the final ejected mass is lost. The last three numbers 
list the final mean metallicities of the cold ISM, the stars 
and the ejected gas. 

The most important result from this table is that the 
final parameters of the models are remarkably insensitive to 
variations of the input parameters. For models DW, the final 
stellar disc mass varies between ~ 4x 10^ and 7x 10^ Af© and 
the gas ejection fraction varies from 0.59 to 0.82. For models 
MW, the final stellar disc mass varies between ~ 2.3 x 10^" 
and 2.8 x 10^" Af© and the gas ejection fraction varies from 
0.12 to 0.22. Figure |l^ shows the evolution of the radial den- 
sity profiles of models MWl and DWl and figure |ll] shows 
the time evolution of the star formation rates, gas fractions 
and gas velocity dispersions. The models of Table 2 behave 
in similar ways, and so these two figures are representative 
of the behaviour of all of the models. These figures are qual- 
itatively similar to those of the simple model of Section 4 
(figure 1^ and ^). The main differences are: 

(a) The gas discs have a sharper outer edge. This is a conse- 
quence of the infall model; the outer edge is determined by 
the final time of the model which sets the maximum cooling 
radius within the halo (c/. figures 

(b) The radial profiles of models MW show oscillatory be- 
haviour near their centres, and the star formation rates and 
gas fractions show oscillatory behaviour as a function of 
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Table 2: Feedback Efficiency: Variation of Input Parameters 



MnHol 


MWl 




DWl 


MW2 


DW2 


MW3 


DW3 


fh 
(Pk. 




0.1 




0.01 




0.1 




^cQ 


0.01 




0.01 




0.03 




^SNII 




0.0 




0.0 




1.0 




^BE 




0.0 




0.0 




0.0 




§5.4 




no 




no 




no 




Ms (Mq) 


2.8 X 10^0 


4.2 X 10^ 


2.3 X lO^*^ 


7.1 X lO'^ 


2.3 X 10^" 3, 


.6 X 10^ 


Mg (Mq) 


5.4 X lO'J 




1.3 X 10« 


5.4 X 10^ 


9.9 X lO'^ 


6.9 X 10^ 1 


.3 X 10** 


{Mq) 


4.6 X 10" 




2.6 X 10** 


7.9 X 10^ 


3.1 X 10* 


6.5 X lO'' 3 


.0 X 10** 


fej 


0.12 




0.59 


0.22 


0.64 


0.18 


0.64 


Tej (Gyr) 


0.25 




0.82 


0.40 


1.8 


0.30 


1.24 


/* 


0.74 




0.10 


0.63 


0.15 


0.63 


0.08 




0.65 




0.03 


0.57 


0.04 


0.56 


0.02 


{Zs/Zq) 


0.55 




0.20 


0.43 


0.30 


0.42 


0.19 


{Z^j/Zq) 


0.29 




0.11 


0.38 


0.14 


0.26 


0.09 



Model 


MW4 


DW4 


MW5 


DW5 


MW6 


DW6 




0.1 




0.1 




0.1 






0.01 




0.01 




0.01 




<5SNII 


1.0 




1.0 




1.0 




<:be 


0.0 




0.05 




0.05 




§5.4 


yes 




no 




yes 




Ms {Mq) 


2.5 X 10^0 


4.3 X 10^ 


2.4 X 10"' 


4.3 X lO'^ 


2.5 X 10^0 4, 


,5 X 10^ 


Mg (Mq) 


2.1 X 10'' 


3.5 X 10^ 


4.9 X 10^ 


9.5 X lO'^ 


1.9 X 10^ 2, 


,9 X 10^ 


M,j {Mq) 


7.1 X 10^ 


3.4 X 10* 


6.8 X 10^ 


3.1 X 10* 


7.2 X 10^ 3, 


,4 X 10* 


fej 


0.21 


0.82 


0.19 


0.69 


0.21 


0.82 


Tej (Gyr) 


0.30 


1.19 


0.30 


1.19 


0.21 


1.12 


/* 


0.73 


0.10 


0.67 


0.09 


0.74 


0.11 


{Zg/ZQ) 


0.62 


0.03 


0.59 


0.02 


0.64 


0.02 


{Zs/Zq) 


0.47 


0.17 


0.44 


0.21 


0.48 


0.19 


(Ze.j/ZQ)) 


0.26 


0.09 


0.27 


0.10 


0.27 


0.10 



time. Both of these effects are a consequence of the galactic 
fountain. 

In these models, the star formation rate begins to rise 
as the gas disc builds up from infalling gas. As the star 
formation rate rises, the cold ISM is converted efflciently 
into a hot phase and this is either driven out of the halo 
or becomes part of the galactic fountain. In models MW, 
most of the gas that escapes from the system is lost within 
this early ( ^ 0.2 Gyr) period of star formation when the 
net star formation rate is close to its peak of ~ lOM0/yr. 
After about 0.2 Gyr, the temperature of the hot phase in 
models MW settles to ~ 10® K very nearly independent of 
radius (c/ figure ^), and so the galactic fountain cycles on a 
characteristic time-scale of ~ 4 x 10* yr. Models DW behave 
in much the same way as the simpler models of Section 4.3, 
except that infall, by construction, extends over a longer 
period of time. In these models, the escape criterion for the 
wind is satistfied over most of the lifetime of the disc and 
hence the model of a galactic fountain is unimportant. 

We discuss briefly the effects of varying the input pa- 
rameters: 

(j)k- The evaporation rate A/ed has a weak dependence on the 
evaporation parameter (pk (oc , equation ^) and obvi- 
ously decreases as (jik is reduced. However, the temperature 
of the hot component is proportional to (^^"'^^ and hence 
rises as <j>k is reduced. The net effect is that the mass of gas 



ejected is relatively insensitive to (jtk, but the mass of the 
final stellar disc increases as (l>k is reduced. 

Ec©: Increasing this parameter reduces the star formation 
rate in the self-regulating model for a fixed gas surface den- 
sity and velocity dispersion (equations |l5| and ^s]). However, 
a lower past star formation rate leads to a higher gas sur- 
face density which increases the star formation rate. These 
effects tend to cancel and so the models are insensitive to 
variations in tc©. 

esNii: Setting this parameter to unity increases the temper- 
ature of the hot component slightly and hence increases the 
efficiency of feedback. As explained in Section 5.3, energy in- 
jection by Type II supernovae at large vertical scale heights 
will always be small compared to the internal energy of the 
hot phase. 

Ebb: Values of ebe ~ 0.05 have little effect on the evolution. 
Provided that ebe is not too large (so that it does not dom- 
inate the net star formation rate), pressure enhanced star 
formation is self-limiting because it increases the velocity 
dispersion of the cold clouds (reducing the star formation 
rate in the self-regulating model) and converts cold gas to 
hot gas. Both of these effects tend to reduce the net star 
formation rate. 

Pressure response of cold cloud radii: Allowing the cold gas 
radii to respond to the pressure of the hot phase provides 



© 0000 RAS, MNRAS 000, 000-000 



16 G. Efstathiou 
1000 



100 



o 
a 

o 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — l—q 

Model MW 




I. . I 



J I I L_ 



6 8 



1000 



100 Ell 



o 
ft 

o 



S 10 b 



0.1 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — l—q 

Model DW 



W 



I 



_l I I I L_ 



100 



o 



43 



N 

I 

O 



10 !- 



0.1 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 




_l I I I I L_ 







4 



6 8 



100 







10 



15 20 




Figure 10. The left hand panels show the evolution of the stellar (dashed lines) and cold gas (solid lines) profiles in models MWl and 
DWl. The right hand panels show the temperature (solid lines) and density (dashed lines) of the hot component. The radius ro is the 
fiducial radius listed in Table 1. The results are plotted at ages of 0.1, 1, 3, 6, 10 Gyr and (for MWl only) at 15 Gyr. 



strong positive feedback in the very early stages of galaxy 
formation when the pressure of the hot phase is high. How- 
ever, most of the cold ISM is ejected on a much longer 
timescale (c/. values of Tej in Table 2) when the typical pres- 
sure of the ISM is similar to that in our own Galaxy. The 
pressure response of the cold cloud radii therefore has little 
effect on the final feedback efficiency. 

The models described here involve a complex set of cou- 
pled equations and a number of parameters. However, one of 
the most interesting aspects of this study is that the equa- 
tions interact in such a way that the evolution of the models 
is insensitive to the parameters. Most importantly, the effi- 
ciency of feedback is insensitive to the thermal conduction 



parameter (p^,. The possible severe suppression of thermal 
conduction by tangled magnetic fields in astrophysical envi- 
ronments is a long standing theoretical problem. However, 
our results show that even a reduction of /t by a factor of 100 
or more will not significantly alter the efficiency of feedback. 

6.2 Chemical evolution 

In this section, we summarize some of the results relating to 
chemical evolution in these models. Our intention is not to 
present a detailed model of chemical evolution in disc sys- 
tems along the lines of, for example, Lacey and Fall (1983, 
1985) but to investigate some of the general features of chem- 
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Figure 11. Evolution of the star formation rate, gas fraction and gas cloud velocity dispersion for the models shown in figure 10. 
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Figure 12. The distribution of stellar metallicities at four radii in model MWl at an age of 15 Gyr. 
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Figure 13. The distribution of stellar metallicities at four radii in model DWl at an age of 10 Gyr. 
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Figure 14. Metallicity gradients in the stars and gas at the final times in models MWl and DWl 



leal evolution with physically motivated models of inflow and 
outflow. The chemical evolution model is based on the in- 
stantaneous recycling approximation as described in Section 
5.6. This is probably a reasonable approximation since the 
timescales of star formation and outflow are ~ 1 Gyr, but 
will overestimate the gas metallicities where the gas density 
is low. As in the previous Section, results from models MWl 
and DWl are used to illustrate the general features of the 
models. The other models listed in Table 2 behave in very 
similar ways. 



6.2.1 Stellar metallicity distribution 

The final mean stellar metallicities are typically Zs/Zq « 
0.5 in models MW and ^ 0.2 in models DW. Models DW 
have a lower stellar metallicity because a larger fraction of 
the ISM is expelled in a wind. The stellar metallicity dis- 
tributions are shown in Figures ^ and (|l|). Figure ^) 
is particularly interesting because this radius is close to the 
solar radius. This metallicity distribution is quite similar to 
that of G-dwarfs in the solar cylinder (see e.g. figure 8.19 
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6.2.3 Effective yields 





isi 

Nl 




According to the simple closed box model of chemical evolu- 
tion, the metallicity of the ISM is related to the gas fraction 
according to 



Zg^-p ln{MJ{Mg + Ms)). 



(50) 



Mg/(Mg + Mj 



Figure 15. The effective yield for model DWl. For each radial 
ring in the galaxy we plot the gas metallicty Zg against the gas 
fraction in that ring. According to the simple closed box model of 
chemical evolution, Zg = —pln{Mg/{Mg + Ms)), where p is the 
yield. The dashed line shows this relation, but using an effective 
yield Peff = p/5. 



of Pagel 1997), showing that the infall model solves the 'G- 
dwarf ' problem of closed box models of chemical evolution. 
The metallicity distributions of model DWl plotted in figure 
(^) also show a lack of stars with low metallicities. 



6.2.2 Metallicity gradients 

Over most of the stellar disc, model MWl has a fairly weak 
stellar metallicity gradient (figure |l^) except at the very 
outer edge where the stellar density and metallicity fall 
abrubtly. This differs from the metallicity gradients seen 
in large disc systems which show linear gradients (see Vila- 
Costas 1998 for a recent review). It is possible that this prob- 
lem might be resolved by including radial gas flows in the 
models (Lacey and FaU 1985, Pitts and Tayler 1985). The 
stellar metallicity gradients in model DWl are steeper, in 
qualitative with observations which indicate that the abun- 
dance gradients in Scd and Irr galaxies are steeper than 
those in earlier type galaxies. 

The radial gas metallicity profiles are shown in figure 
(p^). Model DW contains a gasesous disc extending well 
beyond the edge of the stellar disc. This gas disc has a 
low metallicity in the outer parts, with Z/Zq ^ 10~^ at 
r ^ 2 kpc. At these large radii, the star formation rate is 
always low and the gas disc can survive for much longer 
than a Hubble time without converting into stars. This is 
unlikely to happen in all galaxies for at least two reasons: 
(i) the energy injection from supernovae into this gas will 
not be uniform as assumed in this paper; (ii) the extended 
gas disc is susceptible to external disturbances and so could 
be tidally stripped or transported towards the centre of the 
system in a tidal interaction. Nevertheless, it is possible that 
dwarf galaxies at high redshift possess extended gaseous 
discs, some of which survive to the present day. 



It is well known that the yields derived from applying this 
relation to gas rich galaxies (usually dwarf systems) result 
in "effective yields" , Pcff , that are much lower than the yield 
expected from a standard Salpeter-like IMF. For example, 
Vila-Costas and Edmunds (1992) find effective yields in the 
range p^g ~ 0.004-0.02 and that the effective yield decreases 
with increasing radius. 

The solid line in figure (^5|) shows the final gas metal- 
licity in radial rings in model DWl plotted against the gas 
fraction within each ring. The dashed line shows equation 
( ^o|) with an effective yield of 0.004 {i.e. one-fifth of the true 
yield). The strong outfiows in this model suppress the ef- 
fective yield well below the true yield and produce a strong 
radial variation of the effective yield, in qualitative agree- 
ment with observations. 



6.2.4 Metallicity of ejected gas 

The last line in Table 2 lists the mean metallicity of the gas 
that escapes from the galaxy. The mean metallicity of the 
ejected gas is about 0.3^0 for model MWl and about O.IZq 
for model DWl. In model DWl this value is about 3 to 5 
times higher than the mean metallicity of the final gas disc. 
The ejected gas in this model is therefore 'metal enhanced' 
relative to the gaseous disc. The mechanism for this metal 
enhancement is physically different to that in the models 
of Vader (1986, 1987) and Mac Low and Ferrara (1999). In 
the models of these authors, metal enhancement arises from 
incomplete local mixing between the supernovae ejecta and 
the ISM. In our models, the gas is assumed to be well mixed 
locally, but metal enhancement arises because the gas is lost 
preferentially from the central part of the galaxy, which has 
a higher metallicity than the gas in the outer parts of the 
system. 



6.3 Connection with damped Lyman alpha 
systems 

The column density threshold for the identification of 
damped Lya systems is N{m) > 2 x 10^° cm"^ (Wolfe 
1995) corresponding to a neutral gas surface mass density of 
~ 1.6 Mq/pc^. Comparison with Figure 10 shows that the 
extended cold gasesous discs around dwarf galaxies would 
be detectable as damped Lya systems. Furthermore, in 
CDM-like models, such extended discs around dwarf galax- 
ies would dominate the cross-section for the identification of 
damped Lya systems at high redshift because the space den- 
sity of haloes with low circular speeds is high (Kauffman and 
Chariot 1994, Mo and Miralda-Escude 1994). If this is the 
case, then the metallicities of damped Lya systems would be 
expected to be low at high redshift, Z/Zq ~ few x lO"'^, with 
occasional lines-of-sight intersecting the central regions of 
galaxies where the metallicity rises to Z/Zq ^ 0.1. At lower 
redshifts, the metallicities of damped systems would be ex- 
pected to show a similarly large scatter, but with perhaps 
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Figure 16. The solid line shows the retained baryonic fraction 
1 — /ej as a function of the halo circular speed . The parameters 
^c© eic. adopted are the same as those of models MW4 and 
DW4 listed in Table 2. The dotted line shows the relation adopted 
by Cole et al. (1994), equation (51) with Vf^^t = 140 kms~^ and 
'^hot = 5.5. The dashed line shows equation (51) with v^at = 
75 kms^^ and a^ot = 2.5. 

a trend for the mean metallicity to increase as disc systems 
with higher circular speeds form and the extended gaseous 
discs around dwarfs are disrupted by tidal encounters. 

This is qualitatively in accord with what is observed 
(Pettini et al. 1997, Pettini et al. 1999, Pettini 1999). 
These authors find that the typical metallicity of a damped 
Lya system at 2; ~ 2-3 is about 0.08^0 with a spread of 
about two orders of magnitude. Comparing the metallicities 
of high redshift systems with those of 10 damped Lya sys- 
tems with redshifts z — 0.4-1.5, Pettini et al. (1999) find no 
evidence for evolution of the column density weighted metal- 
licity. Whether these and other properties of the damped 
Lya systems can be reproduced with the feedback model 
described here requires more detailed 'semi-analytic' calcu- 
lations along the lines described by Kauffmann (1996). How- 
ever, the key point that we wish to emphasise is that accord- 
ing to the models described here, most of the cross-section at 
any given redshift will be dominated by largely unprocessed 
gas in the outer parts of galaxies that does not participate 
in the star formation process. The metallicity distributions 
and the evolution of flm as a function of redshift are there- 
fore more likely to tell us about feedback processes and the 
outer parts of dwarf galaxies than about the history of star 
formation. Attempts to use the properties of damped Lya 
systems to constrain the cosmic star formation history {e.g. 
Pei, Fall and Hauser 1999) should therefore be viewed with 
caution. 

6.4 Feedback efficiency as a function of halo 

circular speed and semi-analytic models of 
galaxy formation 

In this section we investigate the efficiency of feedback as a 
function of halo circular speed. We have adopted the param- 



eters of models 4 in Table 1 and run a series of models vary- 
ing the halo circular speed Vv. The virial radius of the halo 
is set to r„ — 150(u„/126 kms~^)^ kpc, the concentration 
parameter c = 10 and the ratio of gas to halo mass within 
the virial radius is set to 1/10. The halo rotation speed is 
set by equation ( |3E| ) with the fiducial disc scale length equal 
to r„/50. With these parameters, the family of models has 
a constant value for the halo spin parameter of Xh = 0.065. 

The retained baryonic fraction, 1 — f^j, is plotted as a 
function of halo circular speed in Figure 16. The dotted line 
shows the relation used by Cole et al. (1994, hereafter C94) 
in their semi-analytic models, 

l-.hot = —^y /?(«„)= f^)""'*", (51) 

1+P(V^) \VhotJ 

where fhot is the fraction of the cooled gas that is reheated 
and Vhot and ahot are parameters. C94 adopt a severe feed- 
back prescription with ahot ~ 5.5 and v^ot ~ 140 kms~^ to 
reproduce the flat faint end slope of the B-band galaxy lu- 
minosity function in a critical density CDM model. The C94 
feedback model does not agree at all well with the models 
described here. There is a slight ambiguity in the appropri- 
ate value of Vv to use in equation (ml) because C94 adopt 
an isothermal rather than an NFW halo profile; the halo 
circular speed at ~ 0.1r„ may be 20% higher than the cir- 
cular speed at the virial radius, but this is far too small a 
difference to reconcile the C94 feedback prescription with 
the models of this paper. 

In fact, the dashed line in Figure 16 shows that our 
models are reasonably well described by equation ( |5l| ) with 
Vhot ~ 75 kms~^ and OLhot = 2.5. Our results therefore sug- 
gest a much gentler feedback prescription than assumed in 
C94. Note that with the C94 parameters, a Milky Way type 
galaxy with Vv ~ 130 kms"'^ would have lost about 60% of 
its baryonic mass in a wind. This is well outside the range 
found from our models for plausible choices of the input 
parameters (c/. Table 2). 

Recently Baugh et al. (1999) and Cole et al. (1999) 
describe semi-analytic models applied to A-dominated CDM 
cosmologies that employ a gentler feedback model. The pre- 
scription for their reference model is similar to that of equa- 
tion ( ^ ) with Vhot = 150 kms"'^ and ahot = 2.0, but with v^ 
replaced by the disc circular speed Udisc- This model is closer 
to the results found here. Assuming angular momentum con- 
servation, a halo with Ah ~ 0.06 will produce a disc with a 
circular speed Wdisc ~ 1.7i7i, (c/. Table 1) and so their model 
can be approximated by equation ( |5l| ) with Vhot ~ 90 kms~^ 
and Ohot ~ 2.0. With these parameters, their model gives 
somewhat stronger feedback than found in our models, but 
is well within the range of physical uncertainties. Kauffmann 
et al. (1993) and Kauffmann, Guiderdoni and White (1994) 
also adopt a much less severe feedback prescription than 
C94 in their semi-analytic models. For a detailed analysis of 
the effects of varying the feedback prescription (and other 
parameters) in semi-analytic models see Somerville and Pri- 
mack (1999). 

The change from an Einstein-de Sitter CDM cosmology 
in C94 to a A-dominated CDM model in Cole et al. (1999) 
partly explains why the revised models provide a reasonable 
match to observations using less efficient feedback. However, 
the revised models predict a faint end slope for the B-band 
luminosity function that is consistent with the observations 



© 0000 RAS, MNRAS 000, 000-000 



Supemovae Feedback in Galaxy Formation 21 



of Zucca et al. (1997) but not with those of other authors 
{e.g. Loveday et al. 1992, Maddox et al. 1998). (The ear- 
lier paper of Cole et al. 1994 attempted to reproduce the 
flat faint end slope of the Loveday et al. luminosity func- 
tion) . The observational differences in estimates of the faint 
end slope of the optical luminosity function are not properly 
understood and so it remains unclear whether a gentle feed- 
back model, of the type proposed here and used in Cole et 
al. (1999), can account for galaxy formation in CDM-type 
models. 

Table 3: Dependence of Feedback Efficiency 
of Model DW on Halo Angular Momentum 



fcoll 


Vc (km/s) 


\h 
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With the Cole et al. (1999) parameterization the effi- 
ciency of feedback depends, by construction, on the surface 
density of the galaxy and hence on the angular momentum 
of the parent halo. In their model, higher angular momen- 
tum haloes lead to more efficient feedback because they form 
low surface density discs with low disc circular speeds. This 
is not what is found in our models. Table 3 lists the ejected 
gas fraction as a function of the halo spin parameter \h. 
Here, the hak) circular speed and virial radius, Vv and r„, 
are the same as for model DW in Table 1, but the amplitude 
of the halo rotation speed (or equivalently the parameter 
fcoll) is adjusted to change the spin parameter of the halo. 
The parameters of the feedback model are the same eis those 
for model DWl in Table 2. The feedback efficiency depends 
weakly (and non-monatonically) on Xh, and is greater in 
systems with low values of Xh- This is because higher sur- 
face densities in low A galaxies result in higher star forma- 
tion rates and a higher temperature hot component that can 
escape more easily from the halo. 

The timescale for feedback in C94 and Cole et al. 1999 
is closely linked to the star formation timescale which is as- 
sumed to be shorter in galaxies with high circular speeds. 
This is not what is found in the models of Table 2. The 
timescale for star formation is of order several Gyr in models 
MW (which have a roughly constant star formation rate at 
late times, see Figure 11), yet the ejection of hot gas occurs 
only in the initial stages of formation with a characteristic 
timescale of ~ 0.3 Gyr. In models DW, the situation is re- 
versed with star formation occuring on a somewhat shorter 
timescale than that for outflow. 



7 CONCLUSIONS 

The main aim of this paper has been to show that super- 
novae driven feedback can operate in a quiescent mode and 
that high rates of star formation are not required to drive 
efficient feedback. In dwarf galaxies feedback occurs on an 
infalling timescale and so can extend over a period of ~ 1 
Gyr. In the feedback model developed here, cold gas clouds 
are steadily evaporated in expanding supernovae remnants 
and converted into a hot component. Critically, the rate at 



which cold gas is evaporated can exceed the rate at which 
mass is converted into stars. If the temperature of the hot 
component is high enough, a wind will form and the hot 
gas can escape from the halo (provided the interaction with 
infalling gas can be ignored) . If the temperature of the hot 
component is not high enough for it to escape from the halo, 
it will cool and fall back down to the disc in a galactic foun- 
tain. Some characteristic features of the models are as fol- 
lows: 

(i) In a Milky Way type system, feedback from supernovae 
may drive out some of the gas from the halo in the early 
phases of evolution {t ^ 0.3 Gyr) when the star formation 
rate is high and the temperature of the hot phase exceeds 
~ 5 X 10^ K. For plausible sets of parameters, perhaps 20 - 
30% of the final stellar mass might escape from the galaxy. 
At later times, the temperature of the hot phase drops to 
T ~ 10*" K and the evaporated gas cycles within the halo in 
a galactic fountain. 

(ii) In a dwarf galaxy with a circular speed --^ 50 kms~^, 
expanding supernovae remnants can convert the cold inter- 
stellar medium efficiently into a hot component with a chac- 
teristic temperature of a few times 10^ K. This evaporated 
gas can escape from the halo in a cool wind. Typically, only 
about 10% of the baryonic material forms stars. Gas ac- 
creted from the halo at ^ 1 Gyr forms an extended gaseous 
disc which, according to the self-regulating star formation 
model used here, can survive for longer than a Hubble time 
without converting into stars. 

(iii) The feedback model developed here is meant to provide 
a sketch of how feedback might operate in a multi-phase in- 
terstellar medium. The model contains a number of obvious 
over-simplifications. For example, we have neglected any in- 
teraction of the outflowing gas with the infalling medium, 
we have not addressed the origin of the cold cloud spectrum, 
ignored the dense molecular cloud component of the ISM 
and neglected any local dissipation of supernovae energy in 
star forming regions. These effects, and other processes, are 
undoubtedly important in determining the efficiency of feed- 
back. Nevertheless, the simplified model presented here con- 
tains some interesting features. Firstly, the model shows how 
positive feedback (via pressure induced star formation) and 
negative feedback (via outflowing gas) can occur simultane- 
ously. Secondly, the models are remarkably insensitive to un- 
certain physical parameters, in particular, thermal conduc- 
tion would need to be suppressed relative to its ideal value 
by factors of much more than 100 to qualitatively change 
the model. If thermal conduction is highly suppressed, it 
may be possible to construct a qualitatively similar model 
to the one presented here in which cold gas is converted into 
hot gas in shocks. 

(iv) The self-regulated star formation and feedback models 
described here provide physically based models for the star 
formation timescale and feedback efficiency as a function of 
the parameters of the halo. The star formation timescale 
and feedback efficiency (or timescale) are taken as free pa- 
rameters in semi-analytic models of galaxy formation (e.g. 
Cole et al. 1999, Kauffmann et al. 1994) and are critically 
important in determining some of the key predicted proper- 
ties of these models, for example, the faint end slope of the 
galaxy luminosity function and the star formation history 
at high redshifts (see e.g. Somerville and Primack 1999). It 
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is therefore important that we develop a theoretical under- 
standing of these parameters (as attempted crudely here) 
and also that ways are found to constrain these parameters 
observationally. The results presented here show that super- 
novae feedback is much loss effective than assumed in some 
earlier semi-analytic models (Cole et al. 1994, Baugh et al. 
1996) but is closer to the more gentler feedback prescriptions 
used in more recent models (Cole et al. 1999, Somerville and 
Primack 1999). 

The feedback model described in this paper has a num- 
ber of consequences and raises some problems which are 
summarized below. 

(i) Evidence for outflows: According to the models described 
here, outflows with speeds of ~ 200 T^g^kms"'^ should be 
common in high redshift galaxies. There is evidence for an 
outflow of ^ 200 kms^^, a mass loss rate of ~ GOMo/yr 
and a star formation rate of ~ 4QMq /yr in the gravitation- 
ally lensed Lyman break galaxy MS1512-cB58 (Pettini et 
al. 2000). The outflow velocity in this galaxy is consistent 
with our models, but the star formation and mass loss rates 
(which arc highly uncertain) arc high. The most likely ex- 
planations are either that MS1512-cB58 is a massive galaxy 
driving an outflow that will remain bound to the system, 
or that it is a less massive system undergoing a burst of 
star formation. In addition to direct detection of outflow- 
ing gas, winds may have other observational consequences. 
The winds from dwarf galaxies will cool rapidly (see Ap- 
pendix B). Wang (1995b) has suggested that photoionized 
gas clouds formed in the cooling wind might contribute to 
the Lya forest. Nulsen, Barcon and Fabian (1998) suggest 
that outflows caused by bursts of star formation in dwarf 
galaxies might even produce damped Lya systems. 

(ii) Damped Lya systems: The extended gaseous discs that 
form around dwarf galaxies in our models have low metal- 
licities because they have low rates of star formation. If this 
is correct, then this largely unprocessed gas would dominate 
the cross section for the formation of damped Lya absorbers. 
The metallicities of most of these systems would be low, but 
would show a large scatter because some lines of sight will 
pass close to the central regions of galaxies containing gas 
of high metallicity. This is broadly in agreement with what 
is observed. Extended gaseous discs would be vulnerable in 
tidal interactions. Some of the gas might be stripped and 
some might be transported into the central regions to be 
converted into stars and hot gas. The evolution of fini de- 
termined from damped Lya systems {e.g. Storrie-Lombardi, 
McMahon and Irwin 1996) might have more to do with in- 
fall, feedback and tidal disruption than with the cosmic star 
formation history. 

(iii) Angular momentum conservation: In hydrodynamic 
simulations, gas is found to cool effectively in sub-units dur- 
ing the formation of a protogalaxy. These sub-units lose their 
orbital angular momentum to the halo as they spiral towards 
the centre and merge. Hence the gas does not conserve an- 
gular momentum during the formation of a massive galaxy. 
(Navarro and Benz 1991, Navarro and Steinmetz 1997, Weil 
Eke and Efstathiou 1998, Navarro and Steinmetz 2000). In 
fact, in the absence of feedback it has proved impossible to 
form discs with angular momenta similar to those of real disc 
galaxies starting from CDM initial conditions. In the models 
described here, it has been assumed for simplicity that the 



specific angular momentum of the gas is conserved during 
collapse. This assumption could easily be relaxed. However, 
the feedback model decribed here suggests that the numer- 
ical simulations miss some important physics. Firstly, it is 
preferentially the low angular momentum gas, infalling in 
the early stages of evolution, that is most likely to be ejected 
from a developing protogalaxy. Secondly, supernovas driven 
feedback may help to solve the angular momentum momen- 
tum problem by ejecting gas efficiently from sub-units. The 
ejected gas may then infall at later times when the halo 
is less sub-structured, approximately conserving its angu- 
lar momentum (Weil et aZ.1998, Eke, Efstathiou and Wright 
2000). 

(iv) Implementing feedback in numerical simulations: There 
have been a rmmber of attempts to implement supernovae 
feedback in gas dynamical numerical simulations {e.g. Katz 
1992, Navarro and White 1993, Navarro and Steinmetz 
2000). These involve either heating the gas around star 
forming regions (which is ineffective because the energy is 
quickly radiated away) or reversing the flow of infalling gas. 
An implementation of the feedback model described here 
is well beyond the capabilities of present numerical codes. 
It would require modelling several gas phases, a cold inter- 
stellar medium, a hot outflowing medium and an infalling 
component, including mass transfer between each phase. It 
might be worth attempting simpler simulations in which cold 
high density gas is added to the halo beyond the virial rar 
dius at a rate that is determined by the local star formation 
rate. 

(vi) Starbursts vs quiescent feedback: It is likely that star- 
burts are more common at high redshift because of the in- 
creased frequency of galaxy interactions. Starbursts could 
contribute to supernovae driven feedback in addition to the 
quiescent mode described here. However, at any one time, 
our models suggest that the cold gas component will have 
a mass of only 20 - 50% of the the mass of the stellar disk. 
Even if a substantial fraction of this gas is transported to- 
wards the centre of a galaxy in a tidal encounter (see e.g. 
Barnes and Hernquist 1996) and is subsequently ejected in a 
superwind, this mode of feedback will be inefficient because 
the mass of gas involved is a small fraction of the total gas 
mass ejected in the quiescent feedback mode over the life- 
time of the galaxy. 

(vii) Metallicity ejection: The mean metallicity of the gas 
ejected from a dwarf galaxy is typically about Zq/W in our 
models, and comparable to the mean metallicity of the stars 
in the final galaxy. Yet typically a dwarf galaxy is predicted 
to expel 5 to 10 times its residual mass in stars. Dwarf galax- 
ies can therefore pollute the IGM with metals to a much 
higher level than might be inferred from their stellar con- 
tent. The high metallicity of gas in the central regions of 
clusters ~ Zq/3 {e.g. Mushotsky and Loewenstein 1997) 
may require a top-heavy IMF and gas ejection from massive 
galaxies. 
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APPENDIX A: FAST COMPUTATION OF THE 
ROTATION CURVE OF A THIN DISC 

We begin with the expression for the potential of a thin disc 

at 2; = 

4){r,z = G) = -2TiG I I Joikr)Joikr')r'fi{r')dr'dk.{Al) 
Jo Jo 

The integral over fc is a well-known discontinuous integral 
(e.g. Watson 1944) 

T tl.\r (l- 'Wt 2 r {l/r)K{r'/r) r' < r 

where K is the complete elliptic integral of the first kind. 
Differentiating equation (A2), we find 



, Ejr'/r) 
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/< (r, r')r' jj.{r')dr' 



J>(r, r')r' fi{r')dr' 



I>{r,r') = 



K{r/r') r'E{r/r') 



r{r'' - r2) ' 



(A3a) 

(A3b) 
(A3c) 



where E is the complete elliptic integral of the second kind. 
This integral is convergent in the limit e ^ 0. The functions 
7< and 7> can be evaluated once and stored, reducing the 
computation of i;^(r) to a simple integral over the surface 
density of the disc multiplied by the tabulated functions. We 
evaluate the epicyclic frequency by differentiating equation 
(A3) numerically. 



APPENDIX B: STEADY SPHERICAL WINDS 

The equations governing a steady spherically symmetric 
wind are 



J.2 



pvr 



dv 

2 2p 



dp d$ , . 

+ pv— = n-c, 

dr 



(Bla) 
(Bib) 
(Blc) 



where q{r) is the mass density injected per unit time and 
TC and C are the heating and cooling rates per unit volume 
{e.g. Burke 1968, Ifolzer and Axford 1970). We assume that 
the gravitational force is given by the NFW halo potential 
(equation 3), d^/dr = vjj{r)/r, and rewrite these equations 
as two dimensionless first order equations 



df_ _ 1_ 

dx 2nx'^{c'^ 



d^_ _ -1 

dx 67^x^(0^ — v^)v^ 
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where c is the adiabatic sound speed, x — r/rv, and all 
velocities are expressed in units of «„. The quanities 7 and 
K in these equations are related to the mass injection and 
cooling rates according to 



j{x,v) 



q{r)M{r) 



k{x, c) = 



Mir)A{T)nl 



M{r) ^A-kI q{r)r dr 



(B2c) 



and the injected gas is assumed to have a uniform initial 
isothermal sound speed of Ci = (fcTi/(0.61mp))^''2. 

We illustrate the behaviour of the wind solutions by 
studying two regimes. Firstly, we assume that q = beyond 
an initial radius Vi — 0.04r^ defining the base of the flow {i.e. 
two disc scale lengths for fcoii = 50). Equations (B2) do not 
have a transonic point when q = (Wang 1995a, see also the 
discussion below) and so we begin the integrations at a Mach 
number slightly greater than unity with — 5Ci/3. We 
adopt the parameters of models MW and DW given in Table 
1 and integrate the equations (Bl) adopting A/ = IOMq/jt 
for model MW and M = 0.2Mo/yr for model DW. These 
mass injection rates are close to the maximum rates at times 
t ~ Tej for the models described in Section 6. The curves in 
Figure 



Bl 



show solutions for initial isothermal sound speeds 
of 0.75, 1.0 and 1.25 times the escape velocity from the centre 
of the halo {vesc ~ 430 kms~^ for model MW and 107 kms~^ 
for model DW). 

The figure shows that the criterion « V2.5ci ^ Vesc 
is about right if the wind is to reach beyond the virial ra- 
dius. For Ci « Vesc the wind in model MW begins at a high 
temperature of Ti ~ 1.4 x 10^ K and cools almost adiabati- 
cally initially, reaching a temperature of ~ 1.5 x 10^ K at the 
virial radius. The timescale for the flow to reach the virial 
radius, ~ 2 x lO^yrs, is slightly longer than the cooling time 
at r„. The behaviour of models DW is quite different. For 
d ~ Vesc the initial temperature of the gas is Ti « 8 x 10^ K 
and cools to ^ 10 K by r = 0.3ri,. As expected from the 
discussion in Section 6, cooling is important in outflows from 
dwarf galaxies (see e.g. Kahn 1981, Wang 1995a, b). 

An investigation of transonic solutions of equations (B2) 
require a model for q{r). An example is illustrated in Figure 
p33 for model DW, using 



M{<x) 

'3' = ^ — 3-exp -r/r»). 



0.04r„. 



(B3) 



In this solution, M(cx)) = O.2M0/yr and the central gas 
density was adjusted to obtain a critical solution for the 
case Ci = Vesc - If the gas is to escape from a dwarf galaxy the 
transonic point must occur before cooling sets in. For such 
systems, the wind parameters would adjust so that a sonic 
point exists at a characteristic cooling scale height as shown 
in Figure B2 . The wind will then cool radiatively just beyond 
the sonic point forming a cold wind as discussed above. It is 
also likely that the wind will be heated to a temperature of 
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Figure Bl. Steady wind solutions for models MW and DW including radiative cooling. The curves show the wind velocity (figs a,c) 
and adiabatic sound speed (figs b, d) assuming that the flow begins at = 0.04r„ with a Mach number of unity. The numbers give the 
initial isothermal sound speed in units of the escape speed tigsc from the centre of the halo. These curves are for a mass injection rate of 
lOM0/yr for model MW and O.2Af0/yr for model DW. 



T ~ 10*K by photoionizing radiation from the galaxy and 
the general UV background. These sources of heating have 
not been included in the models of Figures 



Bl 



and 



The wind will be thermally unstable when cooling sets 
in, and may form clouds. However, in the absence of a con- 
fining medium, the clouds would have a filling factor of order 
unity so the wind is likely to maintain its integrity until it 



meets the surrounding IGM. The external pressure required 
to balance the ram pressure of the wind is 



Pext 



M 



k \ 0.2MQ/yr J \ W kpc 



lOOkm/s 



cm""* K, 



(B4) 
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which is about equal to the pressure of the IGM with a 
temperature of IC^K and an overdensity of 

\l + zj * \^O.2M0/yry 

f ^ \ ^ f '"^ 
\ 10 kpc J \ lOOkm/s 

Provided that the halo is devoid of high pressure gas, the 
cool wind will propagate beyond the virial radius and will 
be halted either by the ram pressure of infalling gas or after 
sweeping up a few times its own mass. As pointed out by 
Babul and Rccs (1992), if a dwarf galaxy is embedded in 
a group or cluster of galaxies with a pressure exceeding ~ 
100 cm~"^K, the bulk motion of the outflowing gas would 
be thermalized in a shock and the cooled shocked gas could 
fall back onto the galaxy generating a new burst of star 
formation. The efficiency of feedback is therefore likely to 
be a function of local environment. 
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